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We have made a search for 57 Fe gamma rays of energy 14.4 keV induced by resonant absorption 
of monochromatic solar axions, as proposed by Moriyama. The proposed axions are suggested to be 
emitted from the Sun, in Ml transitions between the first, thermally excited state and the ground 
state of 57 Fe. An upper limit on hadronic axion mass of 745 eV is obtained at the 95% confidence 
level, it being assumed that z = 0.56 and S = 0.5. 
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The existence of axions, proposed as very light neu- 
tral pseudoscalar bosons which couple weakly to stable 
matter, would solve the strong CP problem associated 
with the 9-vacuum structure of QCD. The effects of 0- 
vacuum lead to an effective interaction which violates CP 
invariance in strong interactions, the magnitude of the in- 
teraction being proportional to the parameter 0. A limit 
on of less than (or of an order of) 1CP 9 follows from 
present experimental bounds upon the electric dipole mo- 
ment of the neutron. In an attempt to guarantee strong 
CP invariance automatically, Peccei and Quinn have in- 
troduced an additional global chiral symmetry, sponta- 
neously broken at an energy scale f a , yielding 0=0 dy- 
namically 0. This solution predicts the existence of a 
pseudo Nambu-Goldstone boson, called the axion, whose 
mass (m a ) is inversely proportional to f a Axions also 
arise in supersymmetric and superstring theories and are 
candidates for dark matter in the universe. In the stan- 
dard axion model [0J^] the Peccei-Quinn (PQ) symmetry 
breaking scale is assumed to be equal to the scale of elec- 
troweak symmetry breaking. The axion mass would be 
roughly of the order of 100 keV to 1 MeV and has been 
experimentally excluded ||. Variant axion models keep 
f a ~ 250 GeV, but drop the constraints of tree-level 
flavor conservation and predict an axion mass near to 
1.7 MeV Q. The existence of variant axions has simi- 
larly been ruled out by experiment ||. In efforts to re- 
tain the PQ idea, new axion models have been proposed 
which decouple the scale of PQ symmetry breaking from 
the electroweak scale and introduce f a at a value much 
greater than 250 GeV. Because coupling constants of ax- 
ions with matter and radiation are inversely proportional 
to f a , axion models of this type are generically referred 
to as invisible axions. The two classes of invisible ax- 
ion models which have been proposed are referred to as 
KSVZ or hadronic axions § and DFSZ or "GUT" ax- 
ions [Q. While astrophysical and cosmological consid- 
erations constrain the mass of these invisible axions to 
a rather narrow range of 10~ 5 eV < m a < 10~ 2 eV, 
but with large uncertainties on either side, less model- 
dependent laboratory experiments have excluded masses 



greater than about 10 keV |p|-|10||. However, it should 
be noted that for hadronic axions there exists a small 
window 10 eV < m a < 20 eV, between the supernova 
(SN) 1987A cooling and axion burst arguments |PHil|, 
as long as the axion-photon coupling is sufficiently small 
(E/Nfn2). 

Because of axion coupling to the nucleus, Moriyama 
|l2| has proposed the existence of almost monochro- 
matic axions produced in the solar interior during Ml 
transitions between the first, thermally excited state of 
14.4 keV and the ground state of 57 Fe . The stable iso- 
tope of iron, 57 Fe (with natural abundance 2.2%), is ex- 
ceptionally abundant among the heavy elements in the 
Sun (solar abundance by mass fraction 2.7xl0 -5 ). The 
axions arc Doppler broadened due to thermal motion of 
the axion emitter in the Sun and therefore they are able 
to excite the same nuclide. In the laboratory it is possi- 
ble that this effect could be measured by absorption of 
solar axions in a 57 Fe target. The possibility of decay of 
the axion into two photons during their traversal from 
the Sun to the Earth is insignificant. The resonant exci- 
tation of 57 Fe by axions would be accompanied by sub- 
sequent decays of excited nuclei, either through emission 
of a gamma ray with an energy of 14.4 keV, or through 
emission of an internal conversion electron. By detecting 
these decay products one could make conclusions about 
the mass of the axion. This makes possible an experi- 
mental test of the hadronic axion window, independent 
of the axion-photon coupling. Experiments based on the 
axion-photon coupling have been proposed and results 
were reported by several authors [fl3l. 

Following the calculations in Ref. |12|1 , one can express 
the total rate of excitation per unit mass of 57 Fe in the 
laboratory as 



R r *r = 1.65 x 10^ 
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where r a /r 7 represents the branching ratio of the Ml ax- 
ionic transition, relative to the gamma transition |l4| and 
contains the nuclear-structure-dependent terms (3 and r\, 
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as well as the isoscalar and isovector axion-nucleon cou- 
pling constants go and g^. 
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where k a and fc are the momenta of the photon and the 
axion, respectively, a ~ 1/137 is the fine structure con- 
stant, 8 r-j is the E2/M1 mixing ratio, (iq — i ~ 0.38, 
/i 3 — 4.71, /3 = -1.19 and r? = 0.80. The quantities fi 
and us denote the isoscalar and isovector magnetic mo- 
ments, respectively, while go and g% are related to the 
axion mass in the hadronic axion model ]l5[ ] as follows: 
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The mass of the axion is related to /„ by 

y/z 1.3 x 10 7 



in a = 1 eV 



1 + z / a /GeV ' 



The constants F and D are the invariant matrix elements 
of the axial current, determined from semileptonic de- 
cays, S is the flavor-singlet axial- vector matrix element, 
and z — m u /md ~ 0.56 is obtained in a first order cal- 
culation of the quark mass ratio. The naive quark model 
gives S = 3F — D = 0.58. However, the estimation of z 
and S suffer from large uncertainties and ambiguity, and 
are still poorly constrained parameters. The error on z is 
of the magnitude of a second-order correction || , while 
the experimental value of S, extracted from the polarized 
structure function data, ranges from —0.09 to 0.57 Jl6| . 
Combined reanalyses have been carried out, from recent 
data taken using different targets, including NLO per- 
turbative QCD 17 1^]. A value S ~ 0.5 was estimated, 
with F + D = 1.257 and F/D = 0.575. 

We have searched for a peak corresponding to the 
14.4 keV gamma ray of 57 Fe in a single spectrum mea- 
sured by a Si(Li) detector (SLP-10180-P, ORTEC). The 
energy resolution at this photon energy is 235 eV for 
presently used facilities. The total internal conversion co- 
efficient for nuclear deexcitation of the first excited state 
of 57 Fe is 8.56 Jl9|], and thus the corresponding probabil- 
ity for emission of 14.4 keV gamma rays is 7 = 0.105. We 
have obtained an energy calibration using the radioactive 
sources 55 Fe, 57 Co and 241 Am (IAEA) with a diameter 
of about 7 mm. The target of 57 Fe was a Mossbauer 
absorber, 95% enriched (WISSEL), with a diameter of 
10 mm and a thickness of about 53 /im. A disc of nat- 
ural iron with the same dimensions was used for back- 
ground measurement. The mass of 57 Fe in the enriched 
target is M — 31.53 x 10~ 3 g, and a distance between 



the target and beryllium window of the detector of 3 mm 
was used. Since the attenuation length of the 14.4 keV 
gamma ray is 20 /jm in iron, the average escape probabil- 
ity from the target is £ = 0.35. We have estimated a total 
efficiency for detection of 14.4 keV gamma rays using the 
13.9 keV gamma ray emission from the Am source (ac- 
tivity 2.05xl0 3 Bq) and obtained e = (1.6 ±0.1) x 10~ 2 . 
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FIG. 1. (a) The 14.4 keV gamma ray peak from the 57 Co 
source, (b) Energy spectra measured with the enriched 57 Fe 
(run) and with the natural iron (background) targets, accu- 
mulated for time periods of 5.3 x 10 6 s. (c) Net number of 
counts in the region of the 14.4 keV gamma ray peak in an 
effort to detect axionic excitation of 57 Fe. 



The position and width of the peak have been esti- 
mated by measuring the 14.4 keV gamma ray peak from 
the 57 Co Mossbauer source, Fig. |l|(a). Energy spectra 



2 



have been measured with the enriched 57 Fe target, pro- 
viding sensitivity to monochromatic solar axions, and 
with the natural iron target, Fig. |l|(b). The difference 
between these two spectra is shown in Fig. 0(c). Data 
and background were both accumulated for time peri- 
ods At = 61.343 days. We have detected N 7 = 56 ± 201 
gamma ray events in the energy interval of 481 eV. Equa- 
tions allow us to relate the axion mass to the 
measured value of N~, 
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where C = fi(3F - D + 25) /3 + {F + D)(l - z)/(l + z). 
If the value of S is 0.5, as suggested by recent analyses 
we obtain m 4 = (4.5 ± 16.0) x 10 10 eV 4 . To de- 
termine an upper limit to the invisible hadronic axion 
mass, we have multiplied our statistical error by 1.645 
and added in the value of m 4 to obtain m a < 745 eV 
at the 95% confidence level. It can be noted that in Eq. 
(||) we have omitted a phase space factor (k a /k) 3 . The 
effect of this on the measured axion mass in our exper- 
iment is estimated to be less than 2 x 10~ 3 . We have 
also considered self absorption of monochromatic axions 
with m a < 1 keV in the solar interior, in addition to 
the effect of absorption by the Earth (day/night mod- 
ulation). The mean free path of the axions has been 
calculated to be (27r) 1//2 (cr 2 m + o 2 ah Y^ 2 /an, the effects of 



nuclear recoil (1.9 x 10 -6 keV) and the redshift due to 
the gravitation of the Sun (1.5 x 10~ 4 keV) being con- 
sidered negligible. The quantities a em (~ 10~ 3 keV) 
and <T a b (~ 10 -3 keV^ 10~ 5 keV) represent Dopplcr 
broadening of the 14.4 keV line of 57 Fe at the tempera- 
ture of the emitter and the absorber, respectively, while 
n is the average density of 57 Fe atoms in the absorber. 
The integrated cross section for axion resonant absorp- 
tion is a = 7ro- rr a /r 7 , where er = 2.6 x 10~ 18 cm 2 
is the maximum resonant cross section of 7 rays and 
r = 4.7 x 10" 12 keV is the total decay width of the first 
excited state of Fe. From the above, it is indicated that 
axions of mass m a < 1 keV would escape from the Sun 
with insignificant absorption (< 10~ 7 ). The day/night 
modulation is less than 10 -5 . 

We have performed the first measurement of nuclear 
axion emission from the Sun, deriving an upper limit to 
the hadronic axion mass of 745 eV. The SN 1987A bound 
for hadronic axion mass is 20 eV p| . Note however 
that both limits depend upon astrophysical and parti- 
cle physics parameters. Solar axion emission has been 
estimated using temperature and 57 Fe density distribu- 
tions, obtained within the framework of the standard so- 
lar model J2(J. SN 1987A axion emission was estimated, 
identifying absorption probability for a given coupling to 
nucleons (g a N) within a simplified model of the super- 
nova density and temperature profile. However, statis- 
tics dominate the uncertainties associated with the axion 
burst argument, with only 19 neutrinos being observed 
from SN 1987A. Observations of other supernovae will 



be required if the uncertainty is to be reduced. Recent 
measurements of the particle physics parameters, used in 
obtaining the SN 1987A upper mass boundary, allow es- 
timation of expected events in the Kamiokande II study. 
From Fig. 1 in Ref. |ll[] we have scaled the events, ex- 
pected due to absorption of 10.96 MeV axions into the 
0~ state of 16 0, using c 2 ,, where c = (3F - D + 25) /6 
is the dimensionless isoscalar coefficient. From the de- 
rived value of g a N ~ 2 x 10" 6 and the recent value of 
9aN = 5 x 10~ 8 m a /(leV) we have recalculated the up- 
per limit of hadronic axion mass from the axion burst 
argument to be m a < 40 eV. 

We conclude, by noting that improvements in detection 
of 14.4 keV gamma rays from 57 Fe atoms are possible, by 
seeking increase in the value of the factor M ■ £ ■ e in Eq. 
(5) and by further suppressing the background. This new 
method of investigation appears very promising in deter- 
mining the hadronic axion window, being independent of 
supernova models and the uncertainties associated with 
them. 
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